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ABSTRACT 

The BL Lac object lES 0229+200 (z = 0.14) has been detected by HESS during observa- 
tions taking place in 2005-2006. The TeV spectrum, when corrected for the absorption of 
gamma-ray photons through the interaction with the extragalactic background light, is ex- 
tremely hard, even if the most conservative level for the background is considered. The case 
of IBS 0229+200 is very similar to that of lES 1101-232, for which a possible explanation, 
in the framework of the standard one-zone synchrotron-self Compton model, is that the high- 
energy emission is synchrotron-self Compton radiation of electrons distributed as a power law 
with a large value of the minimum energy. In this scenario the hard TeV spectrum is accom- 
panied by a very hard synchrotron continuum below the soft X-ray band. We will show that 
recent Swift observations of lES 0229+200 in the critical UV-X-ray band support this model, 
showing the presence of the expected spectral break and hard continuum between the UV and 
the X-ray bands. 

Key words: radiation mechanisms: non-thermal — 7-rays; theory — 7-rays; observations - 
galaxies: general 



1 INTRODUCTION 

The gamma-ray extragalactic sky at high (> 100 MeV) and very 
high (> 100 GeV) energies is dominated by blazars, associated to 
radio-loud active galactic nuclei with a relativistic jet closely ori- 
ented to the Earth. The resulting relativistic amplification of the 
non-thermal jet emission makes blazars extreme objects, with ap- 
parent luminosities exceeding in the most powerful sources 10^* 
erg/s (e.g. Abdo et al. 2009) and variability timescales as short as 
minutes (e.g. Aharonian et al. 2007a, Albert et al. 2007). The inter- 
est for these sources is driven by the possibility to get interesting 
clues on the acceleration processes of charged particles in relativis- 
tic flows and the possible use of their TeV emission to characterize 
the poorly known extragalactic background light (EBL) in the in- 
frared, optical and UV regions (e.g. Stecker et al. 1992, Stanev & 
Franceschini 1998, Costamante et al. 2004, Mazin & Raue 2007). 

The great majority of blazars detected in the TeV band are 
BL Lac objects, showing weak or even absent emission lines. Their 
Spectral Energy Distribution (SED) is characterized by two broad 
peaks, generally interpreted as due to synchrotron (the low energy 
peak) and inverse Compton (at high energy) emission (an hadronic 
origin for the high energy peak is invoked in the hadronic models, 
e.g. Muecke et al. 2003). Since, as suggested by the weak emission 
lines, the environment surrounding the jet is thought to be poor of 
soft photons, the IC emission is likely dominated by the scatter- 
ing of the synchrotron photons themselves. The synchrotron-self 
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Compton (SSC) scenario is supported by a large number of obser- 
vational clues. However, some recent observations are difficult to 
accomodate in this scheme. In particular, one of the most intrigu- 
ing problems is represented by sources with very hard TeV spectra. 
Indeed, even considering the lowest level of the EBL allowed by 
galaxy counts (e.g. Franceschini et al. 2008), some BL Lacs (most 
notably, lES 1101-232 and lES 0229+200) have spectra (harder 
then Flv) oc i.e. raising in the vF{v) representation), chal- 
lenging the standard SSC scenario, in which the decrease of the 
scattering cross section with energy inevitably leads to rather soft 
SSC spectra at TeV energies. 

Different possibilities to produce hard spectra, overcoming the 
problem of the reduced scattering efficiency, have been discussed. 
Katarzynski et al. (2006, K06 hereafter) showed that if the elec- 
trons (assumed to follow a power law distribution with energy) are 
characterized by a relatively large value of the minimum energy, a 
very hard spectrum can be achieved (f^ oc v^^^). The same hard 
spectrum would be visible in the synchrotron component below the 
soft X-rays. Stecker et al. (2007) pointed out that, under specific 
conditions, shock acceleration could produce power law electron 
distribution much harder than the canonical 7~^. We note, how- 
ever, that even if the electron distribution is hard, the rapid decline 
of the efficiency of the inverse Compton scattering at high energy 
inevitably result in soft spectra. Boettcher et al. (2008) propose to 
overcome the problem related to the decrease of the scattering cross 
section assuming that the TeV emission comes from regions of the 
jet at distances much larger than those considered in standard mod- 
els (< 0.1 pc). At large distances the IC emission should be dom- 
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Date Ti r2 -Bfc F" X^/d-o.f. 

7/8/2008 0.7t? 3 1.93 ±0.15 0.85 ±0.25 1.00 23/19 
8/8/2008 1.8 ± 0.1 - - 1.03 23.6/21 



Table 1. Fit of the XRT data of 7/8/2008 (8/8/2008) with the broken power 
law (power law) model. deabsorbed flux in the 0.3-10 keV band in units 
of 10"^" erg cm~^ s~^. In both cases the hydrogen column density is 
fixed to the Galactic value, Nh = 9.2 X 10^° cm"^. 

inated by the scattering of the extremely soft photons of the CMB 
(e.g. Tavecchio et al. 2000), for which most of the scatterings hap- 
pen in the efficient (Thomson) regime. In this model, the X-ray 
and the TeV emission are decoupled, since the X-rays are assumed 
to originate in the inner "blazar region". Finally, Aharonian et al. 
(2008) propose that the absorption of the TeV continuum by a very 
narrow (black body-like) distribution of optical-UV photons in the 
source could explain the observed hard TeV spectra. However, there 
are no indications for the presence of such a very narrow radiation 
field in the environment of BL Lac objects, though the required 
field intensities are within the available upper limits. 

In this letter we address the problem of the hard TeV spec- 
tra for the specific case of lES 0229-1-200 exploiting new data in 
the UV-X-ray band obtained by Swift. The UV-X-ray energy range 
is particularly critical for the model proposed by K06. Indeed, a 
sharp cut-off is robustly expected in the synchrotron emission at 
these frequencies, associated to the low-energy cut-off of the elec- 
tron energy distribution. We use ho = Q.a ~ 0.7, Q,m = 0.3. 



2 THE SED OF lES 0229-K200 
2.1 Swift data 

Swift (Gehrels et al. 2004) observed lES 0229-1-200 two times, on 
2008 August 7 and 8, with a net exposure of 3.2 and 3.3 ksec, re- 
spectively. XRT (Burrows et al. 2005) and UVOT (Roming et al. 
2005) data have been reduced using the standard procedure (de- 
scribed in, e.g., Ghisellini et al. 2007), with the latest version of the 
HEASOFT and CALDB package^ XRT data have been fitted with 
XSPEC12. The data of the first observation are better fitted with a 
broken power law model with a very hard low energy slope (Table 
[T]l, while the second data set is well fitted with a single power law. 
However, above 1 keV the slope and the flux level are very similar. 
For simplicity, in the SED (Fig. [TJ we only report the data corre- 
sponding to the second observation. UVOT magnitudes have been 
corrected for Galactic absorption using the values of Schlegel et al. 
(1998) for the V, B and U filters and the formulae by Pei (1992) 
for the UV filters and converted into fluxes following Poole et al. 
(2008). 



2.2 TeV data 

The HESS spectrum has been measured in 2005-2006 by Aharo- 
nian et al. (2007b, red circles in Fig. [TJ. For such distant blazar 
detected at such high energies, the correction of the observed TeV 
data for the photon-photon opacity is most relevant. 
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In the case of lES 0229-1-200, its redshift implies that the en- 
ergy photons detected at 10 TeV mainly interact with cosmic back- 
ground photons with wavelength Xmax — 1.24:{E.y[TeV]) /im. 
The maximum observed energy of TeV photons is ~ 1 1 TeV, corre- 
sponding to 12.5 TeV in the source rest frame. So the bulk of the ab- 
sorption is produced by background photons around 13-15 fim. At 
these wavelengths the IR background is not directly measurable be- 
cause of the exceedingly large intensity of the inter-planetary dust 
foreground. However the properties of the cosmic sources of these 
photons are known in great detail from three independent survey 
missions: IRAS (observing at 12 and 25 /im), ISO with maximal 
sensitivity at 15 nm, and Spitzer at 8 and 24 /im. Also various 
independent approaches can be used to estimate the background 
intensities. The simplest and most robust one is to determine the 
source number counts over the widest flux density range at various 
wavelengths and to integrate them to get the local background in- 
tensity. The other solution is to calculate the multi-wavelength lu- 
minosity functions at different cosmic epochs (redshift intervals), 
a quite more complex procedure, but suited for computing the red- 
shift evolution of the background photon density. Note that all these 
procedures yield lower limits to the background intensity. 

We have then considered two background levels for de- 
absorption correction. The first one uses the "Low-SFR" model of 
Kneiske et al. (2004), providing minimal IR background level and 
absorption correction. TeV data corrected in such way are reported 
in Fig. 1 as green triangles. We have also computed a corrected 
spectrum (orange squares) using the observationally-based model 
of Franceschini et al. (2008), which takes into account all avail- 
able data on the background sources. The model includes also some 
corrections related to the evolution of the background intensity as 
mentioned above (the latter are small, however, for our source's 
redshift, and tend in any case to reduce the absorption correction). 

The corresponding intrinsic spectrum, even with the conserva- 
tive EBL level assumed in the models, is very hard (photon index 
~ 1.5 ± 0.2) and the high energy peak in the SED appears to be lo- 
cated above 10 TeV. Particularly critical appears the last data point 
(at energy S ~ 11 TeV): once corrected with the Franceschini et 
al. (2008) model, providing a much larger optical depth at these en- 
ergies than the Kneiske et al. (2004) model, the point lies at a flux 
that it is very difficult to reconcile with the current emission models 
(in a similar way to Mkn 501 around 15-20 TeV, see e.g. Aharonian 
2001). However, the very low statistics at the highest energy bins 
and the uncertainty on the EBL spectrum at 10-20 microns do not 
allow yet any reliable conclusion. Some caveats also apply to the 
spectral shape at smaller energies. Indeed, although the spectrum 
is well represented by a power law, at a closer look it hints at a 
more complex shape (this is already visible in the deabsorbed spec- 
trum shown in Aharonian et al. 2007b). In particular, the first bins 
seems to follow a relatively soft power law, while the datapoints 
above 2 TeV describe a harder spectrum (or even a bump). This 
feature in the average spectrum could be either due to variability in 
the spectral shape or could indicate that more than one component 
contributes to the total emission. However, it also depends on the 
specific shape of the EBL spectrum between I and 10 micron. For 
simplicity, and because this feature is not significant, in the follow- 
ing we assume that the intrinsic source spectrum is represented by 
a single hard power law. 

2.3 The SED 

We assembled the SED of lES 0229-1-200 (Fig.[B using historical 
optical-IR data (from NED and Urry et al. 2000), the simultaneous 
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Figure 1. Spectral energy distribution of lES 0229+200. Grey points are from NED. Tlie grey bow-tie is the BeppoSAX spectrum (Donato et al. 2005). The 
open magenta triangle (Urry et al. 2000) is the HST measure for the host galaxy, while the filled magenta triangle is for the core. The red optical-UV and 
X-ray data are from Swift. The red arrow shows a rough Fermi upper limit derived from Abdo et al. (2009). The red filled circles show the observed HESS 
spectrum (Aharonian et al. 2007b), while the green triangles and the orange squares show the spectrum corrected for absorption by interaction with the EBL 
calculated with the LowSFR model of Kneiske (2004) and the model of Franceschini et al. (2008). The black dotted line shows the emission from the host 
galaxy, approximated with a black body function. The other black lines show the synchrotron and inverse Compton models discussed in the text. For the model 
shown by the solid black line we also report the contribution to the total emission of electrons in small energy ranges (see text for details). 



(de-reddened) optical-UV and X-ray data from Swift and the de- 
absorbed HESS data points (average of measures taken in 2005- 
2006). The red arrow shows the Fermi upper limit inferred from the 
absence of lES 0229-1-200 from the list of high-confidence sources 
detected during August-October 2008, the period of the first three 
months of observations, reported in Abdo et al. (2009). 

As shown in Urry et al. (2000) using //5rdata, the optical flux 
is dominated by the host galaxy (open magenta triangle), while the 
flux of the core lies almost an order of magnitude below that of 
the host (filled triangle). Therefore, it is very likely that also the 
IR-optical data from the 2 MASS survey (open gray squares) are 
dominated by the host galaxy. This possibility seems to be con- 
firmed by the UVOT data (red triangles), showing a rapid decrease 
of the emission from the optical to the UV band. We can there- 
fore conclude that the optical (and partly also the UV) datapoints 
are dominated by the emission from the galaxy (approximated as 
a black body in Fig. [T] dotted line) and that the optical-UV non- 
thermal continuum of the jet lies below the UVOT datapoints. 

Unfortunately the HESS data are not simultaneous to the Swift 



pointing. However, as discussed by Aharonian et al. (2007b), the 
TeV spectrum of lES 0229-1-200 showed a quite limited variability 
during the observations in 2005-2006 (small variability also charac- 
terizes the other TeV blazar with a very hard spectrum, lES 1101- 
232). Therefore we think that we can reasonably assume that the 
TeV spectrum is similar to the 2005-2006 average HESS spectrum. 



3 MODELING THE SED 

We reproduce the overall SED with the one-zone synchrotron SSC 
model fully described in Maraschi & Tavecchio (2003). Briefly, the 
emitting region is a sphere with radius 7? moving with bulk Lorentz 
factor r, with a tangled and uniform magnetic field B. The viewing 
angle of the observer is 6, which, in turn, means a Doppler factor 5. 
In the original model the electrons are assumed to follow a smooth 
broken power law energy distribution with normalization K and 
indices n from 7min to 7brGak and n2 above the break up to 7max- 
For some of the applications in this paper (specifically, models 2 
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and 3 below) we also use a single power law, with slope n and 
limits 7min and 7max and normalization K (extrapolated at 7 = 1). 

A "standard" SSC model with a broken power law electron en- 
ergy distribution extending down to relatively low energies (model 
1 in Table|2j fails in reproducing the hard TeV spectrum. The main 
difficulty is to locate the peak of the SSC bump at energy as large 
as few TeV, as required by the data. As detailed in Tavecchio & 
Ghisellini (2008), if one requires a large separation between the 
synchrotron and the SSC peaks a very large Doppler factor is nec- 
essary. Even with 5 = 50 (see model 1 in Table |2]l, the SSC peak 
lies well below 10 TeV. Note also that, in this case, in order to not 
overproduce the observed UV-optical flux we are forced to use a 
rather hard electron distribution, n — 1.5, harder than the canon- 
ical value of 2 (possibilities to have such hard spectra have been 
discussed, see e.g. Stecker et al. 2007 and Virtanen & Vainio 2005). 

The shape of the SED in the optical-X-ray region shows that 
the synchrotron continuum falls quite rapidly from the X-ray to the 
UV band. As already mentioned this is what we would expect in 
the case of a truncated electron energy distribution, with the emis- 
sion being the low-energy tail of the electrons at the lowest energy. 
It is thus tempting to model the entire SED with such a model. 
The solid line in Fig. [T] shows the model assuming a single power 
law with a large minimum electron Lorentz factor (model 2 in Ta- 
ble |2j, 7min = 8.5 X 10^. As discussed in K06, below the typi- 
cal frequency of the electrons at 7min the synchrotron spectrum is 
very hard, showing the characteristics slope Fj, oc v^^^. The corre- 
sponding SSC spectrum shows a peculiar "triangular" shape, with 
a low energy power law reflecting the F^, oc u^''^ slope of the syn- 
chrotron target photons, and a fast decrease after the maximum due 
to the rapid decline of the scattering efficiency. 

To better visualize the situation, we also report (colored 
curves) the synchrotron and the SSC emission of electrons in small 
"slices" of energy, with Lorentz factors in the range 8.5 — 20 x 10^ 
(magenta), 2 - 4 x lO'' (cyan), 4-8x10'' (red) and 8 - 16 x lO'' 
(green). As can be seen, the SSC emission is almost entirely dom- 
inated by the emission from the low energy end of the distribu- 
tion (magenta), also responsible for the hard tail of the synchrotron 
component. While the luminosity of the synchrotron emission from 
electrons of different energies is almost equal, the decrease of the 
scattering efficiency makes the luminosity of the SSC emission to 
strongly decrease for electrons with increasing energy. 

The position of the SSC peak in the Klein-Nishina regime re- 
flects the energy of the electrons at the low energy end, vic ~ 
^iniii{mc? /h)5 (e.g. Tavecchio et al. 1998). The corresponding 
synchrotron frequency is Vs — 3.7 x 10^ Bj^-^^S. Then the ratio 
between the magnetic field and the Doppler factor: 

1 = 1.7x10-^-^:^^^ G (1) 

where iJa.keV = hus is measured in keV and Eic.tcV ~ huic, the 
energy of the SSC peak, is in TeV. Setting Es ~ 0.1, Eic ~ 10 
and S = 30, we find B — 0.5 mO, much smaller than what is 
usually found in other TeV BL Lacs (i.e. B ~ 0.1 — 1; e.g. Celotti 
& Ghisellini 2008). As detailed in Tavecchio & Ghisellini (2008), 
to increase the separation between the two peaks (i.e. move the 
SSC peak at larger energies and/or the synchrotron peak at lower 
energies) we have to increase the Doppler factor. As an example of 
this, we report with the short dashed line another possible model 
with the same SSC peak frequency and a lower synchrotron break 
frequency. In this case the required value of Doppler factor is 5 = 
50. 



4 DISCUSSION 

The combination of the optical-UV and X-ray data of lES 
0229+200 suggests a rapid increase of the flux between the UV and 
the X-ray band, that we interpret (following K06) as due to a large 
value of the high-energy cut-off in the electron energy distribution, 
leading to a very hard synchrotron continuum, oc u^^^. This 
synchrotron spectrum translates into a similarly hard SSC spec- 
trum, reproducing the intrinsic (de-absorbed) TeV data. We note 
that, at a closer look, the spectrum predicted by the model seems to 
be harder than what required by the data derived using both the con- 
sidered models for the EBL. A better agreement would be attained 
by considering larger levels for the EBL, such as that adopted in 
the "Best model" of Kneiske et al. (2004). However, the detection 
of blazars at large redshift at TeV energies (Aharonian et al. 2006, 
Albert et al. 2008) seems to exclude such high levels for the EBL. 

We remark that, within the framework of the one-zone SSC 
model, the solution discussed here is the only one able to reproduce 
hard TeV spectra, without invoking unrealistically large Doppler 
factors. As we have shown, this solution has also a direct conse- 
quence on the synchrotron continuum, predicting a sharp roll-off 
of the emission below the X-ray band. The UV-X-ray data from 
Swift, though not simultaneous with the TeV spectrum of HESS, 
strongly support our scenario. 

The model has some important consequences for the accelera- 
tion process and the role of the energy losses on the electrons. This 
is related to the requirement of an of the electron energy distribu- 
tion with a sharp low-energy end. If there were electrons below this 
limit, even with a very hard distribution, the resulting SSC spec- 
trum would be completely different: the peak, dominated by the 
emission of the low-energy electrons (see Tavecchio et al. 1998), 
would shift at energies well below 10 TeV (see also K06). There- 
fore it is essential to assume that (i) the acceleration process is 
able to push almost all the electrons at energies above 7min and 
(ii) that the electrons, once injected into the emission region, do 
not cool rapidly. In fact, the cooling would inevitably lead to the 
formation of the standard 7^^ power law below 7min. The require- 
ment (i) is particularly demanding for the models based on rela- 
tivistic shocks. This could indicate that electrons are powered by 
reconnection of magnetic fields (e.g. Giannios, Uzdenky & Begel- 
man 2009 for an application to blazars). This conclusion is sup- 
ported by the extremely low value of the magnetic field derived 
in the emission region. Assumption (ii) is instead naturally ac- 
complished in our model: with the parameters used for model 2 
we derive a (observed) cooling time of the electrons at 7mm of 
tcooi 2i 5 X 10**/7miiiB^(5 ~ 10* s ~ 47 months. Therefore, such 
a configuration could be maintained for some years. This cooling 
timescale is valid if only radiative cooling are efficient. Realisti- 
cally, the emitting source is expanding, thus adiabatic losses are 
more efficient in cooling the electrons. The stability of the TeV 
flux for long timescales suggest than that the acceleration/injection 
mechanism can act for long times at a constant level. 

A direct consequence of this scenario is the very low level 
of the emission in the GeV band. Thus a prediction is that Fermi 
will not detect lES 0229-F200. It is possible (or even likely) that 
there are other active regions in the jet, possibly producing 7-rays 
in the Fermi range. In this case, since the emission in the TeV and 
in the GeV band would not come from the same location, a strict 
correlation between the variations in the two bands is not expected. 

Our model for the emission of lES 0229-1-200 reproduces both 
the hard UV-soft X-ray continuum and the very hard TeV spectrum. 
In the already mentioned model of Boettcher et al. (2008), instead. 
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7min 7b 7max « "2 B [G] K [cm R [cm] (5 

1 10'' 6 X 10^ 3 X 10^ 1.5 3.4 8.5 x lO'^ 6 10^'^ 50 

2 8.5 x 10^ - 4 X 10^ 2.85 - 5 x lO-* 3.5 x lO'-* 5.4 x 10^"^ 30 

3 5 X 10^ - 4 X 10'^ 2.85 - 4 x 10""' 6.7 x 10* 5.4 x 10^'^ 50 



Table 2. Input parameters for the models shown in Fig[T](model 1: long dashed line; model 2: solid line; model 3: short dashed line). See text for definitions. 



the UV-X-ray and the TeV components are not related, being emit- 
ted in very different regions of the jet (X-rays in the blazar region, 
TeV emission in the large scale jet). Therefore a very effective way 
to distinguish between these two alternatives is through variability 
of the optical-X-ray and TeV emissions. Related variations would 
point to a common region for the two emissions, challenging the 
model of Boettcher et al. (2008). On the contrary, totally indepen- 
dent variability would indicate that two regions are active, ruling 
out our scenario (and any one-zone model). Moreover, Boettcher et 
al. (2008) predict a larger GeV flux than predicted in our scenario, 
coming from the SSC emission of the inner region responsible for 
the X-ray continuum (albeit the GeV emission could be accom- 
modated with a different choice of parameters). Another possible 
criticism could be that in order to fit the data, 5 = 25 is required 
in the external regions of the jet, contrary to the VLBI indications, 
suggesting very small apparent velocity at parsec scales (e.g. Finer, 
Pants & Edwards 2008). The same problem could afflict our model: 
a strong deceleration from F — 15 — 30 in the emission region 
to r ~a few at parsec scale would imply a radiative luminosity 
comparable to the jet power (see also Ghisellini, Tavecchio & Chi- 
aberge 2005). On the contrary, the derived jet power is more than 
one order of magnitude larger than the radiative power. Unfortu- 
nately there are no measurements of the jet speed in lES 0229-1-200 
at VLBI scale. 

Concluding, the model presented here is capable of reproduc- 
ing successfully the main features of the observed SED, even if im- 
plying extreme physical parameters, in particular for the minimum 
electron energy and the magnetic field. These conditions are rather 
different than those commonly derived in other TeV BL Lacs. This 
could indicate physical conditions in the jet more complex than 
what assumed in the simple one-zone model adopted here, requir- 
ing more sophisticated models (multi-zone, multi-population of 
particles). 

A rather interesting point concerns the consequences of the 
observed TeV spectrum for the modeling of the EBL. Even using 
the strict lower limits coming from the galaxy counts the derived in- 
trinsic spectrum is exceptionally hard, with the point at the largest 
energy being the most problematic. Further observations of such 
hard spectra, in this and other sources, particularly if extended at 
large energies (> 10 TeV), could then be particularly relevant, pos- 
sibly leading to the exciting evidence that something fundamental 
is missing in our understanding of the emission physics of blazars 
or of the propagation of 7-ray photons in the Universe (e.g. De An- 
gelis et al. 2009). 
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